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Abstract. The continued fraction expansion ol the real number x — ao + Xo, 
a G 2, is given by < x n < 1, x' 1 = a n+1 + x„+i, a n+1 G IN, for n > 0. 
The Brjuno function is then B(x) = J^^lo x ° Xl ■ ■ ■ x n-i ln(x^ ), and the number 
x satisfies the Brjuno diophantine condition whenever B(x) is bounded. Invariant 
circles under a complex rotation persist when the map is analytically perturbed, if 
and only if the rotation number satisfies the Brjuno condition, and the same holds 
for invariant circles in the semi-standard and standard maps cases. In this lecture, 
we will review some properties of the Brjuno function, and give some generalisations 
related to familiar diophantine conditions. The Brjuno function is highly singular 
and takes value +oo on a dense set including rationals. We present a regularisation 
leading to a complex function holomorphic in the upper half plane. Its imaginary 
part tends to the Brjuno function on the real axis, the real part remaining bounded, 
and we also indicate its transformation under the modular group. 



1 Hamiltonian Chaos and the Standard Map 

The simplest non trivial model for Hamiltonian chaos is a two dimensional 
real map, called the "Standard Map". It has been introduced more or less 
independently by Chirikov and Taylor [[l]j|] . The occurrence of chaos was dis- 
cussed by Greene || who displayed many numerical results on this model, 
which describes a simplified version of the non-linear coupling of two oscil- 
lators. It occurs naturally in many domains of physics, including celestial 
mechanics, classical quasiperiodic systems, quantum quasicrystals, adiabatic 
response in non-linear mechanics, magnetic toroidal configurations in plasma 
physics, non-linear electronic devices, and many others. 

The Standard Map is a map from the cylinder TxEto itself, defined as 

K K \ 

+ r + — sin(27r6») (mod 1) , r + — sin(27r6») . (1) 

Z7T 2tT I 



Note that the second variable can also be taken modulo 1, in which case we 
a map T 2 -> T 2 . 1 

Hamiltonian form 

2tt 



get a map T 2 — > T 2 . The above map can be written in two equivalent forms 



r' -r = — sin(27T0) , 0' - 6 = r' , (2) 
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• Lagrangian form, where one considers now the twice iterated map, that 
is (0, r) -> (9\ r') -» (9", r"), so that 

9" -20' + 9=— sin(27r<9') , (3) 

Z7T 

and for the n-th iterated map one gets 

9 n+l - 29„ + 9 n ^ = — sin(27T0 n ) . (4) 

Z7T 

This last equation is sometimes called the Frenkel-Kontorova model |Q which 
describes equilibrium positions of a chain of material points placed in a peri- 
odic potential and submitted to an harmonic elastic force between two neigh- 
bouring points. 

If K = 0, we get the so-called "twist-map", which gives 

r n +i = r n = r = p = constant , n +i = 9 + np (mod 1) , (5) 

after n iterations. This map is nothing but a rotation of angle 2irp — we say 
that the rotation number is p. The orbits are all transverse to the axis of the 
cylinder. They are made of a finite number of points (and therefore discrete) 
if p is rational, and they are dense on transverse circles if p is irrational. The 
cylinder is sliced along orbits at irrational values of r, which are intertwined 
with discrete orbits at rational values of r. The question is : in which sense 
is such a pattern stable under perturbations, that is here when K ^ 0? 

Among the orbits which are dense in a curve wrapped around the cylinder, 
particularly interesting are the orbits which will persist under perturbation, 
in particular because they separate the space into domains which do not 
communicate. It is known that when K is large (for example K > 4/3, see 
||), such orbits do not exist, and on the other hand, when K is small, some 
of the irrational orbits persist, depending on arithmetical properties of the 
rotation number. For perturbed twist maps, we define the rotation number 
as liirin^oo n~ l 9 n , where {9 n ,r n ) is the n-th iterated map obtained from (1). 

Other kinds of invariant curves may occur, attached to elliptic periodic 
orbits. For instance if K > is sufficiently small, the point (1/2,0) is an 
elliptic fixed point. Due to KAM Theorem (see |(| for a review), there exist 
homotopically trivial invariant curves on the cylinder winding around this 
fixed point, which form the so-called elliptic islands. We shall not consider 
here the problem of such orbits, although there existence is very important 
in connection with ergodic theory. Indeed one expects chaotic behaviour for 
K large, but the persistence of elliptic islands could prevent the map from 
being ergodic. 

2 The Critical Constants 

For the standard map, we consider now the homotopically non-trivial invari- 
ant curves i.e. wrapped around the cylinder. A natural way to look for their 



Real or Complex Brjuno Functions 3 



existence is to replace the angular variable 6 by the new variable <p £ T 

6 = cj) + u{<t>,K,p) , r = p + u((f>,K,p) -u((f>- p,K,p) . (6) 

With the condition 1 + du/d<fi > 0, it would describe a curve around the 
cylinder, on which the map is expressed as </>' = <p + p, when <fi describes T 
for K and p fixed. We say that (6) expresses on the curve the conjugacy of 
the map to a rotation. The existence of a function it(</>, K, p), analytic in the 
variable <p, insures the existence of an analytic invariant curve with rotation 
number p. We are interested to determine the critical constant K c (p) as being 
the largest possible value of K for which such an analytic function u exists. 
Of course, one could consider regularity constraints weaker than analyticity, 
leading to other critical constants. We look for a perturbation expansion of 
the function it, and we follow the notations of Q. From the standard map 
we get from (6) 

K 

u(<f>+p,K,p)-2u(<p,K,p)+u(<p-p,K,p) = —sm(2ir((f> + u((j>,K,p)) . (7) 



2tt 



For k > 1, we have 



u^^+p,p)-2u^{<i ) ,p)W k) {<P-p,p) = -Uin(27T0 + 27™(0,if,p)) 

llT fc-i 

(8) 

where in the right hand side, one keeps only the terms of order k — 1 in the 
expansion on powers of K. We use now the Fourier series expansion on </>, 
that is v,( k \(j),p) = X^ez uf\p)e 2l7rv ^ , and we see that the coefficient of 
e 2 "^ in the left hand side of (8) is 2(cos(2vri/p) - l)vi k) (p). Therefore (8) 
allows a recursive computation of the Fourier coefficients ui^ (p) , and we get 
for u( k \tp,p) expressions as trigonometric polynomials in (p. However terms 
of the kind 2(003(2^^,0) — 1) occur in the denominators along the steps of 
the recursion. Such factors are called "small divisors" , some of them vanish 
when p is rational, and may become arbitrarily small when v becomes large, 
for irrational p. Now let K c (p) be the minimum over (f> of the convergence 
radius of the expansion 

00 

u(q>,K,p)=Y,K k u^(<j),p) . (9) 

k=l 

For p rational, (8) cannot be solved, and we set K c (p) — 0. For irrational 
values of p, Berretti and Gentile || were able to control K c (p) using the 
Brjuno function B{p) which is a number theoretic function which will define 
in the following Section 4. More precisely, there exists C > such that, for 
any irrational p, 

In ((^(p))- 1 ) -2B(p)\ <C. (10) 
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The functions K c (p) and e 2B ( p > both vanish on all rationals, but the previous 
equation shows that the ratio K c (p)/e~ 2B ^ remains uniformly bounded at 
every irrationals. The fact that this ratio remains bounded is in itself amazing, 
but it recalls earlier and now classical results by Yoccoz |J on the linearisation 
of holomorphic maps. We shall see later that we may have even better results 
in the framework of holormorphic maps. 

The determination of the radius of convergence in (9) is not the whole 
story. It is possible that the function u may be analytically continued for real 
values of 4> and for K > K c {p) real. Thus we would have another critical 
constant K c (p) such that we still have real analytic curves for K c {p) < K < 
K c {p) real. The numerical results |l0| seem to indicate that this is indeed the 
case: see jll| for a detailed discussion of this issue both from the numerical 
and the analytical points of view, which also uses results of ]l^ , |T3t . The 
definitive answer is not known to us today, although we are led to expect 
that the function B(p) plays a central role in the determination of K c (p) (see 
also Davie @). 



3 Complex Analytic Maps 

The problems of the critical constant is better understood in the case of the 
complex analytic maps. We have already seen in (10) that the critical constant 
K c (p) of the complexified version of the standard map is controlled by the 
Brjuno function. A simpler example is the "Semi-standard Map" , which is a 
two dimensional complex map on the cylinder, closely related to the standard 
map (1) : to get the semi-standard map, just replace in (1) the sine function 
sin(27T0) by its positive frequency part (1/2*) exp(2iir8). The procedure to get 
analytic invariant curves proceeds in a completely similar way as Equations 
(6) to (9), and it was proven that in this case [|1|[L6), the critical constant 
Kssm(p) defined in a same way as above, fulfils as in (10) 

\]n((K ssm (p))- 1 )-2B(p)\<C . (11) 

The numerical results (especially the figure 16) in ref. |l6| provide more. Not 
only the ratio K ssm (p)) / e~ 2B ( p ^ is bounded on irrationals, but it is extend- 
able to a continuous function on [0, 1], bounded below and above by positive 
constants. This result is amazing if one remember that both K ssm (p)) and 
e -2B(p) van i sn a t all rationals. Therefore the Brjuno function B(p) is a good 
model to represent the singular behavior of In ( ( J\T ssm (p) ) — 1 ) ■ 

The Brjuno function was introduced by Yoccoz Q in the apparently sim- 
pler problem of the linearisation of complex holomorphic maps around their 
fixed points. This is a more than one century old problem (see jlTj for a nice 
review), which we can state as follows. Let f(z) be a holomorphic map such 
that /(0) = 0, f'(0) — e 2l7rp . Is it possible to conjugate the map / to its 
linear part? This means that we look for a function h, holomorphic in a disk 
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of radius R f , such that h(0) = 0, h'(0) = 1, and f(h(z)) = h(ze 2iT P). Note 
that such a function h, if it exists, is unique. In this case, the function / 
is said to admit a Siegel disk of radius Rf. The Siegel disk is a topological 
disk with conformal radius Rf, since it is the image through the normalised 
conformal map h of the disk \z\ < Rf. 

We quote now the classical results on this question |H| . i) If p is rational, 
there is no disk, that is Rf = 0. ii) if p is irrational and satisfies a (strong) 
Liouville condition, we still have Rf — 0. Hi) if p is a diophantine irrational 
(see Section 5 below), then there exists a Siegel disk and Rf > 0, more 
precisely this happens when B(p) is finite, iv) If B(p) = +oo, then that there 
exist functions / such that Rf — 0. Indeed Yoccoz [|| proved the following 
: define R(p) as the smallest radius of the Siegel disks Rf obtained when / 
varies in the compact family of all univalent maps on the unit disk such that 
/(0) = 0, and /'(0) = e 2l7Tp . Then we have 

lln^p))- 1 )-^)! <C. (12) 

Now consider the family of quadratic polynomial Pi{z) = e 2np (z — z 2 ), and 
call i?2(p) the radius of the Siegel disk associated to it. Observe first that, 
through the rescaling z — > e~ 2lnp Rx z, then Pi is transformed in e 2l7rp z — Rz 2 . 
In the rescaled variable, we see that R2 is the maximum value of the constant 
R for which a circle with conformal radius one is persistant. Therefore Riip) 
is the critical constant adapted to the present case, and this leads to the 
analogies between (10), or its equivalent in the real case, and (12). 

Here again, the numerical results (now the figure 6) in ref. |l6[ | bring 
some continuity properties. Not only the ratio R2(p))/e^ B ^ is bounded on 
irrationals, but it is extendable to a continuous function on [0,1], bounded 
below and above by positive constants. Therefore the Brjuno function B{p) 
is again a good model to represent the singular behavior of In ((i?2(/o)) _1 ) • 
In our work Jl8| which started from these observations, we give arguments 
which strongly support the conjecture that the ratio R2(p))/e~ B ^ is not 
only continuous but satifies a Holder continuity condition with exponent 1/2. 
More precisely, the Brjuno function displays the universal singular behaviour 
(up to some Holder- i continuous function) of the critical functions occuring 
in small divisors holomorphic problems in dimension one. 

The use of the Brjuno function was somewhat implicit in the work of 
Buric et al. p9[, where they attempted to find representations of the critical 
constants by what they called modular smoothing. Singular functions of the 
same type occured in MacKay p^] in relation to the Brjuno condition. 

It is nevertheless useful to recall here briefly one of the steps, called 
renormalisation, which plays a special role in Yoccoz 's argument, and will 
appear to be crucial in understanding the fine regularity properties of ratios 
of the type R 2 {p))/e- B{p \ For this purpose, we follow H), and we we con- 
sider first a rotation of angle 2irp, with < p < 1, that is z — > e 2lnp z, 
acting in an open disk of radius R p centered at the origin in the com- 
plex plane. We need an arbitrary point a, such that \a\ = R pi and for 
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simplicity we take the real point a = R p . Let a' = e 2 ™ p R p its image. 
Now, consider the angular sector bounded by the lines Oa and Oa', namely 
A p — { z | < arg z < 2-np , \z\ < R p }, the line Oa' being excluded. 
Consider the orbit made of the successive iterated points z ni n > 1 start- 
ing from zq G A p , and let z q the first of these points which also belong to 
A p . The map zq — * z q is thus the first return map in the sector. We have 

(27r) _1 (arg z q ) = (27r) _1 (arg zq) +qp—l. We now take in the sector the vari- 

_ i 

able u such that its complex conjugate u — z<> , where now u belongs to a disk 
i_ 

of radius Rp . For the values uq and u q , corresponding to zq and z q , we have 
(27r)^ 1 (arg%) = (2 7 r)- 1 (arg u ) - q+ ± (mod 1) = (2vr)- 1 (arg «o) + - p - The 
original map which acted in a disk of radius R p , leads in the new"renormalised 
variable" u, to a rotation with rotation number - , acting in a disk with radius 
i?(!/ p ), such that lni? p = ph\R(i/ p y 

This construction extends to the non linear perturbed case, for example 
P2(p, z) — e 2np (z — z 2 ), with a lot of complications. Suppose that there is a 
Siegel disk for Pi^p, z). In this disk, there are conformal coordinates on which 
the maps is exactly a rotation of angle 2irp, and on these coordinates we apply 
the linear renormalisation. The problem is then to give an interpretation of 
the renormalised coordinates u which we obtain. It appears that there exist a 
holomorphic map in the variable u with rotation number p~ , which admits 
a Siegel disk, with conformal radius R(i/ P ) such that \nR p — plnRn/ p -\ = 0. 
However, this map is not a polynomial with degree 2. This led Yoccoz to 
extend the problem to the compact family of univalent map on the unit 
disk with rotation number p, and he has considered the minimum R{p) of 
the radius of the Siegel disk taken over this family of maps. The result is two 
modifications to the relation lni? p — p\nR^i/ p ^ = obtained in the linear case. 
First due to the minimisation procedure, the best one could get is a positive 
uniform upper bound for this expression instead of zero. Second, there is a 
special difficulty when p goes to zero. In this case the Siegel disk is strongly 
distorted, since there is an other fixed point which tends to zero when p goes 
to zero. The comparison between the linear and the non linear case becomes 
unjustified in this limit. Yoccoz proved that the result is an additional loga- 
rithmic term in the estimate, so that we only get that In R p — p In R(i/ P ) — bap 
is bounded. It is therefore natural to compare the function — \nR(p) (as well 
as — lni?2(p)) to the solution of the equation B(p) — p\nB(l/ p) + Inp = 
which we will see, is nothing else than the Brjuno function. 

4 Continued Fractions and the Brjuno Function 

We first give a somewhat unusual definition of the continued fraction expan- 
sion sometimes called "japanese continued fractions" |2~]}| . Let a be a fixed 
real number such that ^ < a < 1. Then, given the starting number x, the 
coefficients a n and e„ are recursively uniquely defined by the conditions 

x = ao + £02:0, and Vn > 0, x^ 1 = a n +i + e n+ ix n +i , (13) 
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with Vn > 0, a — 1 < e n x n < a. We define the modified integer part [x] a 
and the modified fractional part {x} a as follows, 

[x] a = [x — a + l]i and {x} a = {x — a + l}i + a — 1 , (14) 

where [x]i and {x}\ are the usual integer and fractional parts of x (so that 
< < 1)- With these notations, we can rewrite (13) as 

ao = [x] a , £qxq = {x} a , and, a n+ i = [z" 1 ]^ e n+ ix n+ i = {a;~ 1 } ce ,Vn > . 

(15) 

Therefore the x n are generated by iterating the function A a (x) = \{x 1 } a \, 
that is Vn > 0, x n+ i = A a (x n ) = ({a^ 1 }^ — jx" 1 — A more de- 

tailed description states that the map A a is made of the following branches 

branch k + : A a (x) = — — k for < x < — , (16a) 

v ; x k + a ~ k 

branch k~ : A a (x) = k for — < x < . (16b) 

When ^ < ol < 1, the function A a maps the interval [0, a) to itself, whereas 
when a = i, it maps the interval [0, a] to itself. In both cases, it is convenient 
to set A a (0) — 0, and we get a map which is infinitely differentiable by pieces, 
and the points where it is not differentiable accumulate to 0. Now x and the 
reduced fraction p n /q n admit the following representation 

£o Vn £o 

x = ao-\ =- , — = ao + 



£i ' a ~ , £l 

ai H In ai + 



£n— 1 • £n— 1 

a 2 + • H a 2 + • H 

(17) 

As long as the x n 's do not vanish, we have 

X = ; , X n = (-£„) , (18) 

q n +qn-l£nX n Pn-1 ~ Xq n -1 

and the recursion relations 

Pn = a n p n -i + e n -iq n -2 , Pa = a , P-i = 1 , (19a) 
q n = a n q n -i + e n -iq n - 2 , qa = 1 , q~i = , (19b) 

so that we get < go < qi < 92 < • ■ • < q n < q n +i < ■ ■ ■■ We also define 

f3 n = X X! ■■■X n = (-l) n £o£l ■ ■ ■ S n (q n X - p n ) , (20) 

and we have 

: < (3 n q n+ l = < - • (21) 



+ a q n+1 + e n+1 q n x n+1 a 
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Now there exist A(a), with < A(a) < 1, and positive constants G\ and C%, 
such that@ 

0n < C\X(a) n , q n > C 2 X(a)- n . (22) 

Indeed we have 

for V ^~ 1 < a < 1 , A(a) = A(l) = = 0.618... , (23a) 

and for i < a < 1 , A(a) = A = >/2- 1 = 0.414... . (23b) 

When x n =0 for some n, and a; m 7^ for to < n, then we have x = p n /q n 
which is rational, and we say that the fraction stops at order n (with our 
conventions, we have x m = 0,Vn > to). Conversely, if x is rational, the 
continued fraction expansion stops at some finite order n. For a — 1, we get 
the classical Gauss continued fraction expansion for which all signs e n = +1, 
and for a = 1/2, we have the continued fraction to the nearest integer. Note 
that when a/1, the results of equations (23a) and (23b) are not obvious. 
For details, and in particular for the extension to others values of a, with 
< a< i, see @|2). 

Given a positive real function / on (0, 1), the Brjuno series Bf(x) is the 
sum (which can be infinite) of the series with positive terms 

00 

Bf"\x) = f} n -if{x n ) = f(x )+x f(xi) + .. . + x xi ■ ■■x n -if(x n ) + . ■ ■ , 

n=0 

(24) 

where | < a < 1, and for k > 0, x% defined in Eq. (13) or (15). As mentioned 
above, when x is rational, we have x n = for some n, and we use x m = for 
m > n, The following results are easily obtained from the definitions 

Bf\x) = B ( f a) (x + l), (25a) 



Bf'(x) = xBf 1 I - I + f(x) for < x < a , (25b) 

B { f a \x) = Bf\-x) for < x < 1 - a. (25c) 

In particular, B ^ {x) is an even function. More surprising is the following 
result g§ : in the a = 1 case, for B { f ±} {x) = %(Bf\x) ± Bf\-x)) which 
are the even and odd parts of B^p {x), we have for < x < i, 



Sj +) (.T) = x J Bj +) (I) + ~G(a;) , with (26b) 



(26a) 



G(x) = f(x) + f(l - x) + (1 - x)/ + 2^5 - • (26c 











+ 2^ } 
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In order to prove the previous equations, we use Equations (23a-c) and the 
succession of transformations 

1 1 1 x l-x l l 

—x — » l — x — > > l = ► = l — > ► x , 

l — x I — x l — x x x x 

which provides the requested relations between B(x) and B(—x). 
Now, it is convenient to introduce the following specific notations: 

i) In Bj (x), when a = l, we omit the superscript (a), and when a = |, 

we replace the superscript (a) by e, so that Bf(x) = Bf{x) and Bj(x) = 
B^ 2 ^ [x] respectively. 

ii) We omit the subscript / when f(x) = — ln(x) = ln(x _l ), so that 
B^ a \x) = B ( ™l(x), B{x) = B [ \{x) and B e (x) = B^£\x) respectively. 
We will call B(x) the Brjuno function, which has been mentioned above in 
Sects. 2 and 3. We have 

oo 

B(x) f3 n -xln(x~ 1 )= - In(af )-afo m(xi)+. . .-x xi ■ ■ ■ x n ~! hi(x n )-. . . , 

n=0 

(27) 

where the x n are obtained from (13) using a = 1 (Gaussian case), whereas 
B e (x) is given by the same equation (27) with x n obtained from (13) using 
a — i (continued fraction to the nearest neighbour). Both functions B(x) 
and B e (x) are 1-periodic, and take value +oo for x rational. From (25a), the 
odd part of B(x) is given for < x < \ by B_{x) = ix(ln(x _1 — 1)), which is 
continuous (and even Holder continuous for any exponent a < 1). Moreover, 
B e {x) is even, and it has been proven |D| that the difference B e (x) — B + (x) 
is not only bounded, but continuous, and even Holder continuous for expo- 
nent -i. This refines a more general statement |l^,|22j which says that the 

differences B^ — B{x) are bounded over the irrationals. 

The numerical computation of B(x) and B e {x) is delicate, due to the 
instabilities of the continued fraction expansion. However, it is very easy to 
compute their values when the continued fraction expansion is periodic, that 
is when x is an irrational quadratic number. This applies to noble numbers, 
in which case the x n are constant after a certain order. 

5 The Brjuno Series and Diophantine conditions 

A real number is said to be a Brjuno number if and only if B(x) is finite, 
and we also say that x satisfies the Brjuno diophantine condition. Brjuno 
numbers are irrationals and real numbers satisfying the classical diophantine 
conditions (which we recall below) are Brjuno numbers. In [H, we show that 
for | < a < 1, B^ a \x) is finite if and only if x is a Brjuno number. More 
precisely, the proof says that for any a £ [0, i], the difference \B^ (x)~ B(x)\ 
is bounded over irrational values of x. We also show that for a = 1, the 
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difference \B(x) — X)^o In 1 l n (9n+i)| is bounded over irrational values of 
x, so that we recover the original definition of the Brjuno numbers p3[| : x 
is a Brjuno number if and only if X^^Lo In 1 m (9n,+i) is bounded over the 
irrational. One can see jl^,|2^] that such a definition of the Brjuno numbers 
does not depend of the particular value of a used to compute the q n . 

We now report the usual definition || of the diophantine conditions : 
we say that x is an irrational diophantine number of order r > (and we 
write x G C(r)), if there exists c > such that for any integers p and q, 
such that q > 0, we have \x — p/q\ > cq~ 2 ~ T . Some classical facts need to be 
recalled here j24j. First, for any p and q such that < q < q n +i, we have 
\qx — p\ > \q n x — p n \, where p n /q n is the Gaussian reduced fraction to x. 
Therefore, in order to have x S C(r), it is sufficient to check that for any 
n > 0, \x — p n /q n \ > cq n 2 ~ T . Second, Liouville's classical theorem asserts 
that algebraic numbers of degree n belong to a; € C(n — 2). Moreover Roth's 
theorem shows that all algebraic numbers belong to C(r), for all r > 0. 
Finally, for an arbitrary irrational, and any n > 0, we have (g n <7„+i) _1 < 
\ x ~ Pn/ln\ < Qn 2 - Using (21) for a = 1, we get an equivalent caracterisation 
of the diophantine conditions: x e C(t) if and only if there exists a constant 
c > such that fi n > c /3 n t \ for any n > 0. 

Now we introduce for v > 0, the Brjuno series B^y(x) = B x -» (x) for the 
fonctions f{x) = x~ u (still using a = 1), 

OO OO / ~ \ — oo 

n=0 n— v 7 n— 

(28) 

Using (21) one gets 

oo oo 

Tl = 



The series -B^^x) converges if and only if the series Qn 1 ^ \Qn x — Pn\~" 

converges, that is if the series J2n°=o ( hi 1 ^ 2v \ x ~ {Pn/<ln)\^ v also converges. 
As a consequence, if Br u y(x) < oo, then q^ 1 ^ 2 " \x — {pn/Qn)\~ U is bounded, 
and x £ C(l/V). Conversely, assume r > 0, and a; £ C(r), then we have 

oo 

B W (x)<c-^ g - 1+ -. (30) 

Using bounds in (22), we get the following statement: If x £ C(r), then for 
any i/ such that t < , Bs v \{x) < oo. Therefore, there is a relation between 
the diophantine conditions C(r), and the convergence of the Brjuno series for 
f(x) = x~ v : the set of irrationals x such that B^{x) = B x -»(x) is bounded, 
is contained in C(l/v), and contains C(— e + for any < e < \jv. 

In some sense, the Brjuno conditions is related to the limiting case v = 0, 
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and in particular, x G C(r) for r > implies B{x) < oo, that is x is a 
Brjuno number. Using a more general function /, positive on (0,1), and 
monotoneously decreasing in the vicinity of zero, we can introduce a wide 
family of conditions Bf(x) < +oo. The diophantine conditions obtained will 
be mainly governed by the singular behavior of / around zero. A power 
law behaviour would simulate the usual conditions, whereas a logarithmic 
behaviour would generate a condition similar to the Brjuno condition. Other 
interesting examples would be obtained by taking functions / of the form 
x~ v \ log(x)\^, x~ v \ log(x)| M | log(| logxDI* 7 , and so one. 

6 The Brjuno Operator 

We will introduce now some functional analysis in order to solve Equations 
(25a-c). For fixed \ < a < 1, let us consider the operator T( a ), acting on 
locally Lebesgue integrable functions / on the real line, which verify 



f(x) = f(x + 1) for almost every x € 1R, (31a) 

f(x) = f(—x) for almost every x e (0, 1 — a) . (31b) 

The operator is defined by 

(T (a) /)(x) = xf (-) , if xG(0,a) . (32) 



It is understood that the function Ti a )f is completed outside (0, a) by impos- 
ing on T( a )f t ne same parity and periodicity conditions which are expressed 
for / in the above equations (31a-b). The functional equations (25a-c) can 
then be written in the form 

(l-T^B^^f. (33) 
This suggest to study the operator T( a ) on the Banach spaces 

X a , p = {/ : IR -» R | / verifies (31a-b) , / e L p (0, a)} (34) 
endowed with the norm of L p (0, a), namely 

||/lkp= \f(x)\> dx^) ' P , (35) 

for p e [1, oo]. Note that one could also use L p (0, 1), instead of L p (0, a), and 
that if p < p' one has the obvious inclusion X a _ p i C X a/p . If (1 — Tr a \) is 
invertiblc in the considered space, then (25a-c) have a unique solution for 
Bf*\ provided that the / in the right hand side of (25b) also belongs to the 
space. The invertibility property is given by the following theorem, which 
states in particular that the spectral radius of T( a ) is strictly smaller than 1 . 
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Theorem. Tr a \ is a linear bounded operator from X Q)P into itself for all 
a G [h, 1] and for all p £ [l,oo]. Its spectral radius on X ap is bounded by 
the constant A(a) of Equation (22), and therefore 1 — 21 a ) is invertible. 
For the proof, see |l8fl . We will just observe here that the result is immediate 
in the p = oo case. Indeed, 

P? a) f)(x) = P n -i(x)f(x n ) = /? n _i(z)(/ o A n a ){x) , (36) 

where the map A a is defined above (see (16a-b)). Therefore 

||T ( " Q) /|U i00 < Bup x (A l _ 1 (a:))||/|| ad , < cA(ar- 1 ||/|| Q , p (37) 

and one gets the theorem (for p infinite) by taking the 1 [n— th root of both 
sides. For the other values of p, it is convenient to make use of the mea- 
sure which is invariant under transformation by the map A a , instead of the 
Lebesgue measure. An immediate consequence of the theorem, is that if we 
take f(x) = ln(x), for < x < a, then / € X a}P , for all finite p and therefore 
we also have B^ G X ap for all finite p. 

However, we have a stronger property in the a = | case. Here, we set 
again T e = Tn/2). In this case, the logarithm belongs also to the set X* C 
X a=1 / 2 .p=i, made of even, periodic (with period 1) functions belonging to 
the so-called "BMO-space" . In this space, / has bounded mean oscillation, 
more precisely the following semi-norm ||/||* is bounded, with 

||/||. = <^P]7j- y I/ — /rldas , (38) 

where the mean value of / over / is fj = /, \f — fi\dx, and the sup 
is taken over all possible intervals I with length |/| smaller than one. The 
BMO space has remarkable properties. First it is contained in all L p spaces 
for p finite, and it contains the space, second it is the space adequate to 
describe functions having singular behaviour not worse than logarithmic, but 
around every point in a dense set of the real line, and third, it has remarkable 



properties connected to the harmonic conjugacy transformation 25 18]. 

In |L8]]2(|, we have shown that 1 — T e = I — Tm/2) is invertible in X*, and 
therefore that B e = B^l^ G X*. Since \B a - B e \ is bounded for \ < a < 1, 
we have the unexpected consequence that all B^ for | < a < 1, also have 
the Bounded Mean Oscillation property, although it cannot be shown directly 
through the properties of 1 — Tf a \ . 

The BMO property obtained for the Brjuno function in the real case, was 
one of our motivations to consider the complexification procedure which we 
will describe in the last Section of this paper p7| . 

7 Application to Holder— continuous Functions 

In this Section, we will consider only the case a = \. In this case the map 
A a = Ax/2 is continuous on the interval (0, |]. The functional equation for 



Real or Complex Brjuno Functions 13 



the Brjuno function B e for a = 1/2 is 

[(l-T e )B e ]( 2 ;) = -logx, (39) 

for all x £ (0, 1/2), complemented with the condition that B e is even and 
periodic. We suppose that the right hand side of this equation is pertubcd, 
by an additional term /, which is less singular than the logarithmic function, 
and we want to study the singular properties of the perturbed solution. Since 
the equation is linear, we only need to consider the action on / of T e and 
(1 — T e ) _1 , which we will conveniently recall the Brjuno operator B e . We will 
consider even and periodic functions / which are continuous. It is sufficient 
to know the value of / on [0, 1/2], so we assume / <G C® 1 j 2 y One can check 
that Tf is also continuous provided we set T/(0) = 0. We need now the usual 
Holder's type semi-norms for continuous functions : let / € C® jy 2 p then wc 
define the Holder's 7-norm as 

\fU= sup '^-y , (40) 

o<x< v <i/2 \x-yy 

with < 7 < 1. This is a seminorm since it vanishes on constant functions, 
so that we introduce the norm: 

H/H^AI/I. + SI/U , (41) 

where |/|oo is the L x norm of /, and A and B are positive constants which 
we can choose arbitrarily, provided that they do not vanish. We say that 
/ G C 7 , if / G Cjq 1 y 2 j and |/| 7 is finite. We now have: 

Proposition. T e is a bounded operator in C 7 , for the norm ||/|| 7 , when 
< 7 < 1/2, provided B/A is large enough: if B/A > (2 7 - 2" 7 ')- 1 , the 
norm of T e corresponding to the norm (41) satisfies ||T e || 7 < 2( 27 ~ 1 ) < 1. 
Therefore for < 7 < 1/2, T e is a contraction, and 1 — T e is invertible. 

We need the following Lemma 
Lemma. Let < y < x < 1/2, and define x\ and yi by the following 
conditions 

1 1 

y = , x = , (42) 

n + yi m + xi 

with n > 2 and m > 2, and -1/2 < x 1 < 1/2 and -1/2 < y x < 1/2, then wc 
have 

IN-NI<^- (43) 

Proof of the Lemma. Since y < x, we have n — m > x\ — y\ > —1. therefore 
n > m. Let n — m = p > 0. We have x — y = xy(p + yi — x\). So that we need 
to prove ||xi| — \yi\ \ < \p + yi — x\\. This is obvious when p = 0. We always 
have \ \xi \ — \yi\ \ < 1/2, so the required inequality also holds when p > 2. In 
the remaining case p = 1, we set 77 = sign(yi) and e = signer), and we need 
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to check that ||a;i| — \yi\\ < + — e|a;i||. Still because the left hand side is 
smaller or equal to 1/2, this last inequality is not obvious only when rj = — 1 
and e = +1. It therefore remains to show that ||a;i| — \yi\\ < |1 — \yi\ — |xi||- 
Setting u = 1/2 — \x\\ and v = 1/2 — \yi\, the last inequality is equivalent 
to |1 — v/u\ < |1 + v/u\, which is readily checked since u/v is real and non- 
negative. 

Proof of the Proposition. Let < y < x < 1/2, and x\ and yi as in the 
preceding lemma. We have 

\T e f(x) - TJ(y)\ = \xf(l/x) - yf(l/y)\ = \xf(\ Xl \) - yf(\ yi \)\ (44a) 
<|x- 2/ ||/(|x 1 |)| + |y||/(|x 1 |)-/(| 2/1 |)| (44b) 
< k-2/||/|oo + |y||/| 7 |ki|-|j/i|r (44c) 

<k-^||/|oo + |/| 7 ^=^, (44d) 
where we have used / G C 7 , and the Lemma. Therefore 

M £:% m * \ x - + (r) 7 m 1-27 ^ ^ 

<(l/2) 1 - 7 |/| 00 + (l/2) 1 - 27 |/| 7 , (45b) 

since < y < x < 1/2, and 7 < 1/2. For y = 0, the right hand side can be 
replaced by its first term (l/2) 1_7 |/| 0O , and the above inequality extends to 
the case where y vanishes, so that 

|T e /| 7 < K 7 (f) = 2^ 1 \f\ oc + 2 27 - 1 !./ 1 7 . (46) 

For the norm, we get 

||T e /|| 7 = AITe/^ + BITe/loo (47a) 
< 2 27 - 1 [A|/| 7 + + 2- 7 A)|/U] , (47b) 

<2( 27 - 1 )||/|| 7 , (47c) 

provided 2~ 27 B + 2^ A < B, that is A/B < 2 7 - 2^ 7 which completes the 
proof. The above proposition has two obvious consequences 

• Since C 7 C C 7 whenever 7' < 7, we have 

/ G C 7 and 7 > 1/2 => TJ G C 1/2 (48a) 
/ G C 7 and 7 > 7o , 7o < 1/2 => TJ G C 70 (48b) 

• When 7 < 1/2, T e is a contraction on C 7 . Therefore 1 — T e is invertiblc 
and B e = J2(f(T e ) n = (1 — T e ) _1 preserves C 7 , and we have 

fEC and 0<7<l/2 => B e / G C 7 (49a) 
/ G C 1/2 B e / G C 7 , V7 such that < 7 < 1/2 . (49b) 
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We have reproduced here the proof of [£6| , because it is essentially elementary. 
In fact we have a slightly better result for B e than for T e , as shown in the 
next proposition, which shows that the C 1 ^ 2 property is effectively reached. 
Its proof JH],|2(| is too difficult to be reproduced here. We have 
Proposition. If / e C 7 , and 7 > 1/2, then B e f € C 1 / 2 . 

The Brjuno function B e which we have studied in the previous section 
is nothing else than H e £, where £ is equal to minus the logarithmic fonction 
restricted to [0,1/2]. When made even and periodic, this function is not 
continuous. Suppose that we perturb I by a function with enough regularity 
properties (for example C 1 or C 1 / 2+€ ), the change in B e £ will be continuous 
and even in C 1 / 2 , that is Holder-^ continuous. In this sense, the 'most singular 
part' of the Brjuno function is stable or 'universal', roughly speaking modulo 
Holder-^ continuous contributions. As noticed at the end of Section 5 above, 
we can use either B e = i^ 1 / 2 ) or B = B^ since a similar argument starting 
from (26a-c) shows that their difference is also Holder-^ continuous flis| . 

This provides a frame to understand the properties of the critical con- 
stants K of the Sections 2 and 3 above. We assume that the singularity 
comes from the renormalisation equation (1 — T e )K — £ + / , and not from 
additional singular behaviour coming from / in the right-hand side. If it 
would exist, such an additional singular behaviour would require a further 
physical interpretation. This argument, which is usual in the renormalisation 
analysis of singularities, was one of the motivation for our previous work 
| fl8| . The renormalisation equation allows naturally to conjecture that the 
difference between the Brjuno function B e , and the logarithm of the various 
critical constants (multiplied by a suitable coefficient), is continuous and even 
Holder-^ continuous. As an example, we conjecture that the ratio of e~ B and 
the radius of the Siegel disk of the quadratic polynomial, is an Holder- ^ con- 
tinuous function of the rotation number. These conjectures are in agreement 
with the numerical results Jl6 10 1 . 



8 The Complexification of the Brjuno Function 

We consider here the case a = 1, and we want to associate to the function 
/ in Xi 2) a function @, holomorphic in the upper half plane, such that 
Im^ — ► Bf when z goes to the real axis. Since Re<£ is associated to the 
harmonic conjugate of /, we expect to find better boundedness properties 
when / has the BMO property. We will here describe our procedure, and 
report the results p7[ . 

We associate to / a function F(z) holomorphic in <D\[0, 1], and vanishing 
at infinity, as follows 

F(z) = - [ ^-dx . (50) 

7T J X - Z 

For x real, we have Imf (1 ± ie) = ±f(x) for x £ [0, 1], and Imf (1) = 
for x ^ [0, 1]. We will be particularly interested in the case f(z) = ln(z), in 
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which case we get F(z) = —it 1 Li2(l/z), where L12 is the classical dilogaritm 
function p8| . Now we set 

+N 

${z) = lim ^F(z + n) > ( 51 ) 

—N 

and we get a function <P holomorphic in the upper-half plane H + , periodic 
with a real period equal to one, and such that for x real, ImF(a;±!e) = ±f(x). 
In fact the previous equation defines a pair on functions <P±, respectively 
holomorphic in the upper or lower half plane H-|-, so that the natural frame 
in which our procedure takes place is the frame of complex hyperf unctions, 
which we will not consider here |^7| . 

We consider now the action of T, with (T f){x) — xf{l/x) if < x < 1, 
/ and Tf being complemented using periodicity. Using the above correspon- 
dence / 1—* Tf, a correspondence F *— > TF is induced on holomorphic func- 
tions in C\[0, 1], vanishing at infinity. We get 

(TF)(z) =-zjr (f (~-m\- F(-m) \ + ^ F'(-m) . (52) 

m=l V V Z / / m —l 

In fact, TF is essentially — z Ylm=i ^ — m )' U P ^° an & ^ nc additive 
correction, which could be determined by the vanishing condition at infinity. 

If / is associated to F as above, the solution Bf(x) (for a = 1) of (25a-c) 
is associated to the series 

00 

^(z)^^™, (53) 

7L m=0 

where we use the notation F for S^^'_ 00 F(z + n) understood as the 
symmetric summation (51) to insure convergence. 

It is now interesting to display the link between (53) and the modular 

group GL(2, It). Let g = e GL(2, 7L), which mean a,b,c,d G 2, 

e g = ad—bc= ±1. To g we associate the following group action on functions 
holomorphic on C\[0, 1], that is g 1— > L g F, with 



{L g F)(z) = (a-cz){F(± 



1 ^ /,(' £\\ _!/,<(' !l) . ,r, 4 ) 



Let Ai+ C GL{2,7L) be the multiplicative monoid generated by the unit 

matrix, and the set of matrices ( ? ^ | , for m > 1 integer. The monoid 

V 1 my 

A4+ C GL(2, 7L) can also be defined as the set of matrices including identity 
and the matrices ( ° \ J G GL(2, Z) such that first, d > c > a > 0, and 
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second, d > b > a > 0. In GL(2, Z) there is a unique product decomposition, 
namely Vg € GL(2, Z) there exist a unique set of three matrices k, m and 
h, with g = fcm/i, and m € A4 + , k E Z, where Z is the translation subgroup 

of matrices ^ ^ ^ , n e Z, and h e H, where /i is the order eight sugroup 

of GL(2, TL) made of the matrices ^^j, and , with e = ±1 and 

e' = ±l. 

Now (53) is rewritten as 

= E E (*) ■ (55) 

The double sum over g and /c amounts to a sum over a part of the full modular 
group (here one over eight). The contribution over the seven other possible 
parts would be —Bf(z), ±Bf(z~ 1 ), ±Bf(—z), and £>/(— z^ 1 ). 
We will now summarize the results: 

i) The sums in (55) converge in the open upper half plane as long as / is 
in Li(0, 1) which insures that F is holomorphic in C\[0, 1], and vanishes at 
infinity. 

ii) When / is in L p (0, 1), p finite, then Bf is in the Hardy Hp space. 

iii) If / is such that F has bounded real part, then the same holds for Bf. 

iv) For f(x) = ln(x), and F(z) = — 7r _1 Li 2 (l/z), we get the complexified 
Brjuno function, B, holomorphic in the upper half plane, vanishing at +ioo. 
When z goes to a real number in a non-tangential way, we have the following 
limits when e > goes to zero : the real part Re B(x + ie) has a bounded limit 
for any real x. This limit is continuous at all irrationals and has a decreasing 
jump of n/q at each rational p/q. When a; is a Brjuno number, lmB(x + ie) 
goes to the Brjuno function B(x). 

The limit properties of the complex Brjuno function on the real axis are 
characteristic of functions / having a logarithmic singularity aroud zero. Al- 
though the boundedness of the real part reminds the BMO property of the 
real Brjuno function, it is in fact a stronger property. This is one more re- 
markable feature of this function. We are convinced that the interpretation 
of the properties of the Brjuno function in terms of the modular group is 
promising. On the other hand, we can hope to find an interpretation of the 
complexified version of rotation numbers analogous with the usual interpre- 
tation of complex frequencies in terms of damped oscillations, but this is 
another story. 
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